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By symbolic summation methods based on difference fields we present a general strategy that transforms definite 
multi-sums, e.g., in terms of hypergeometric terms and harmonic sums, to indefinite nested sums and products. 
We succeeded in this task with all our concrete calculations of 2-loop and 3-loop massive single scale Feynman 
diagrams with local operator insertion. 



1. Introduction 

We aim at the simplification of multi-sums to 
expressions in terms of indefinite nested sums 
and products. The calculations are performed in 
the context of the computation of 2- and 3-loop 
massive single scale Feynman diagrams with lo- 
cal operator insertion. These are related to the 
QCD anomalous dimensions and massive opera- 
tor matrix elements. The mathematical expres- 
sions depend on the Mellin variable N, defined 
either at the even or odd positive integers. The 
Feynman diagrams may be mapped to definite 
nested sums over hypergeometric expressions in- 
volving also harmonic sums |1I2) . Given these 
usually huge expressions, the task is to find sim- 
pler representations that can be processed further 
in physics. At 2-loop order all respective calcula- 
tions are finished |3I4I5) and lead to representa- 
tions in terms of harmonic sums only. This is not 
necessarily the case from 3-loop order onwards. 
In particular, during our calculations, we left the 



'Supported by the Austrian Science Fund (FWF) grant 
P20347-N18. 

t Supported in part by SFB-TR-9 and EU TMR network 
HEPTOOLS 

* Supported in part by SFB-TR-9 

§ Supported by the Austrian Science Fund (FWF) grants 
P20162-N18 and P20347-N18. 



solution space of harmonic sums, and only in the 
end, while combining our results, we end up again 
at harmonic sums or their generalized versions [6] . 

Inspired by symbolic summation of hypergeo- 
metric terms [7], in particular by creative tele- 
scoping [5] and finding hypergeometric solutions 
of linear recurrences jS], we developed a general 
summation machinery for indefinite nested sums 
and products. Based on a refined version [10] 
of Karr's summation theory of n£* -difference 
fields pj] we use algorithms from |12I13|10I14| 
[T5] implemented in the Mathematica package 
Sigma [16] that attacks the summation quantifiers 
from the innermost sums to the outermost sums 
and which transforms step by step the sums un- 
der consideration to indefinite nested sums and 
products. Finally, we use the Mathematica pack- 
age HarmonicSums |17j inspired by various ideas 
of [11216118119] and extensions of it [20] in order to 
transform the indefinite nested product-sum ex- 
pressions -whenever possible- to harmonic sums 
or its generalized versions [BJ. In particular, com- 
bining the results, the elimination of algebraic re- 
lations among the occurring harmonic sums plays 
a crucial role. 

For the physical background and further results 
in our computations we refer to [21] . 
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2. Simplification of multisums 

We consider the following problem: Given a 
definite nested multi-sum over hypergeometric 
term^l involving also harmonic sums, for instance 

JV-2 J+l -j+N+r-2 

( _ ir+s (J +i)(-W- 2 )(-, +j V-2)!H 

(JV-«)(*+l)(-j + JV + r)! ' 1 j 

/md an alternative representation in terms of in- 
definite nested product-sum expressions. This 
means in particular indefinite nested sums over 
hypergeometric terms (like binomials, factorials, 
Pochhammer symbols) that may occur as polyno- 
mial expressions in the numerator and denomina- 
tor with the additional constraint that the sum- 
mation index ij of a sum X)/=\ f(ij) mav oc- 
cur only as the upper index of its inner sums 
and products, but not inside of the inner sums 
itself; for a precise but rather technical definition 
we refer to |15j . Typical examples are harmonic 
sums |l|2j defined by 

S mi ,...,m k (N) = 

N - ifc — l 

E sign(mi)' 1 sign(m fc ) Ifc 

• |mi| / i -l^fcl 

»i=l H ifc=l l k 

with mi,...,mfc being nonzero integers, their 
generalized versions called 5-sums [5] 

S mi ,...,m k • ■ ■ i Xki N) — y ] ~fmi ' ' ' ^ ] Int i 

ii=l 1 i h =l %k 

(3) 

with constants Xi or, e.g., nested binomial sums 
of the form 

N ^ 7 N ^ 7 

i=l V i / i=l V i / 

5 /GO is hypergeometric in j iff f(J + l)//(j) = r(jr') for 
some fixed rational function r(j). 



examples for these sums are given in Section [3] 
Note that the sum (Q} is not given in indefinite 
nested form since, e.g., the summation index j 
of the outermost sum occurs inside in the sum- 
mands of the inner sums; for a transformation to 
indefinite nested sums see Subsection 12.21 

2.1. The underlying summation principles 

In order to transform sums such as |T]) to indef- 
inite nested sums, we use the summation pack- 
age Sigma [16] by the following strategy; the 
underlying algorithms |12ll3HOI14ll5j are based 
on various extensions and refinements of Karr's 
summation theory of n£*-difference fields [TTj . 
We process each summation quantifier and trans- 
form step by step the occurring sums to indefinite 
nested versions. Suppose, e.g., that we derived 
already an expression f(N,j) in terms of indefi- 
nite nested sums and products w.r.t. j and that 
we want to attack the next summation quantifier, 

sa y J2fj 2 f( N ,j)- If f( N J) is free 01 N , we are 
already done. Otherwise, we proceed by creative 
telescoping. 

Deriving recurrences by creative telescoping 
Given an integer d > and given a sum 

a 

F(a,N) :=^2f(N,j) (5) 

with an extra parameter N, find constants 
co(N), . . . , Cd(N), free of j, and g{N,j) such that for 
< j < cl the following summand recurrence holds: 

c (N)f(N, + + c d (N)f(N + d,j) 

= g(N,j + l)-g(N,j). (6) 



Creative telescoping has been originally intro- 
duced with Zeilberger's algorithm |8I22) where 
the summand f{N,j) can be a hypergeometric 
term and the derived Ci(N) are rational func- 
tions in N and g(N,j) is a rational multiple of 
f(N,j). More generally, with the summation 
package Sigma the summand f(N,j) can be an 
indefinite nested product-sum expression w.r.t. 
j, and one searches for constants Cj(iV) which 
are indefinite nested product-sum expressions in 
N and expressions g(N, j) which are indefinite 
nested product-sum expressions w.r.t. j. 
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If one succeeds in this task to compute a sum- 
mand recurrence, one gets by telescoping the re- 
currence relation 



c (N)F(a,N) 



-c d {N)F(a,N + d) 
g(N,a + l)-g(N,0). 



(7) 



Remark. In most instances, a depends linearly 
on N, which means that g(N, a + 1) itself might 
be a definite sum. In this case, one first has 
to transform g(N,a + 1) to an indefinite nested 
product-sum representation by recursive applica- 
tion of our strategy. In other instances, a has to 
be sent to infinity to deal with sums of the form 

If one succeeds in this task, we obtain 



co(N)F(N) 



■c d (N)F(N + d) = h(N) (8) 



where the Cj(iV) and h(N) are given in terms of 
indefinite nested product-sums. Then we can ap- 
ply another feature of Sigma. 



Recurrence solving 



Given a recurrence of the form J8|, find all solutions 
in terms of indefinite nested product-sum expressions 
(also called d'Alembertian solution). 

Based on the underlying algorithms, see e.g. [23] 
[53], the derived d'Alembertian solutions of © 
are highly nested: in the worst case the sums 
will reach the nesting depth d — 1. In order 
to simplify these solutions, a refined telescop- 
ing paradigm [15] based on depth-optimal IIS*- 
diffcrence fields [TD] is activated. This machin- 
ery delivers algebraic independent sum represen- 
tations [Tl] with minimal nesting depth [T3] . 

To this end, one tries to combine these solutions 
such that the evaluation agrees with the original 
sum F(N) for the first N = 0, ...,d- 1 initial 
values. Then since both expressions are a solution 
of a recurrence of order d, the expressions agree 
-up to some mild side conditions- for all N > 0. 

As worked out in this subsection, Sigma is a 
general symbolic summation toolbox that assists 
in the task to transform definite multi-sums to 
indefinite nested product-sums. At this point we 
should mention that the strategy sketched above 
might fail: it might go wrong to represent f(N, j) 



in {5| in terms of indefinite nested sums or prod- 
ucts, the Ansatz of creative telescoping might fail 
to find a recurrence of the form ([5]), or we might 
miss sufficiently many solutions (i.e., some of the 
solutions of ([SJ are not expressible in terms of in- 
definite nested product-sum expressions) so that 
they can be combined to an alternative represen- 
tation of the input sum ([5]). However, for all ex- 
amples that we have encountered so far in 2-loop 
and 3-loop computations, the machinery never 
failed. 

For the transformation of our results to 
harmonic sums and their generalized versions, 
J. Ablinger's HarmonicSums package [17 plays an 
important role. It is inspired by various ideas 
of |1I2I6I18I19| and extensions of it [20]; some 
highlights are given in Section [3] 

2.2. A detailed example 

Subsequently, we illustrate this interplay of the 
different Sigma-functions for telescoping, creative 
telescoping and recurrence solving in the differ- 
ence field setting by the given sum ([I]) which can 
be rewritten in the form 

3=0 



-< (-j + N + r)\ 



j+N+r-2 (_iy r-j+N+r-2 



E 

s=0 



(N-s)(s + l) 



(9) 



We start with the innermost sum 

-j+N+r-2 , . \ s (-3 + N+r-2\ 

F (N,j,r)= J \,^ s , J (10) 



s=0 



(JV-«)(s + l) 



and compute the recurrence 

(-j + N + r-l)F(N,j,r)+ 
(j-r + l)F(N,j,r + l) 



-j + N + r 



Next, we solve the recurrence and find the general 
solution 

(-j + N-l) r , 1 



(-i-l)r {N + l)(-j + N + r-l) 

with the generic constant c; here (a)k denotes 
the Pochhammer symbol defined by (a)k = a(a + 
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1) . . . (a + k — 1) for k > 1 and (a)o 
we determine the first initial value 

F o (N,j,0) = - 



1. Finally, 



0'-iv + i)(iv + i) 



(iV - l)AT(iV + l)(2 - A/%- 

by applying again our summation strategy. 
This tells us that we have to choose c = 

(-l) w (j + l)! 
(JV-l)iV(JV+l)(2-JV)j 



so that both sides of 
1 



(JV+l)(-j' + JV + r-l) 
(-l)^(j- + l)!(-j + JV-l) r 
(N — 1)N(N + l)(-j - l) r (2 - AO, 

agree for r = 0. As a consequence, the identity 
holds for all r > for generic values j, AT. By care- 
ful checks it turns out that it holds for all AT, j, r 
with N > 2, < j < N - 2 and < r < j + 1. 
Summarizing, we transformed the definite hyper- 
geometric sum pop to a linear combination of hy- 
pergeometric terms. At this level, we could have 
used also the well-known hypergeometric machin- 
ery presented in [7]. 

Now we turn to the double sum in ([9]) which 
we already simplified to a definite single sum with 
hypergeometric terms involved: 



Fi(N,j) 



E 



=J i-j + N + r) 



F (N,j,r). 



Again, we compute a recurrence with Sigma: 
(-j-2)F 1 (N,j)-F 1 (N,j + l) 



3 



N 



(N + l)(-j + N-2)(N-j)\ 

(-!)"(,• + 2)(j 



1 



D! 



(N - 1)N(2 - N)i 

Solving this recurrence, we end up at the follow- 
ing alternative representation 

F 1 (NJ) = (-iy(j + l)\x 

1 ( BflHj) , jv 2 +i 

N \ \ (N + iy(-j+N-l) (N-1)N(N+1)* 
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(-!)' 



N- 



1 ^ 

i=l 



N-i)\(i+l)\(N- 



(11) 



in terms of indefinite nested product-sums in j. 
Now we are ready to deal with the triple sum © 
itself which is given in the form 

N-2 

F ( N )= 2)1^1 (N,j) 

3=0 

where the summand Fi(N,j) is given by (|11|) . 
Creative telescoping in the general setting of dif- 
ference fields leads to the following recurrence 

(N + l)F(N) + (-N - 2)F(N + 1) = 
A^ 3 - 4N - 4 



N 2 {N + 1)(N + 2) 2 

(-l) 7 ^ 4 + 6N 3 + 5N 2 - 2N - 4) 



(N ~1)N 2 {N + 1) 2 {N + 2) 2 



(—1) N N\ 



N-2 

E 



(AT+ 1)(AT + 2) f-* H+^R(iWRW 

where the right hand side can be simplified fur- 
ther to 

(N+l)F 2 (N) + (-N-2)F 2 (N+l) = + ^)_ j 



-iV 3 - 4N 2 - 8AT - 4 



A-2(Af + 2) 2 
S x (iV) 



iV 2 (AT + l) 2 (N + 2) 2 (A" + 1)(N + 2) 

applying our strategy again. To this end, we solve 
the recurrence and find the closed form 



F(N) 



-N 2 -N-l (-1) N (N 2 + N + 1) 



(12) 



N 2 (N + 1) 3 N 2 {N+l) 3 
Si(N) S 2 (N) 2S_ 2 (AQ 
(N + l) 2 N + l N+l 

of the triple sum (fTJ). 

3. Some recent calculations 



Subsequently, we present some selective exam- 
ples from our recent computations [21]. Here 
several thousand multi-sums have been simpli- 
fied to indefinite nested product-sum expressions 
following the strategy from above. In order to 
carry out these computations mechanically (i.e., 
by pressing just one button), a new package 
called EvaluateMultiSums has been developed 
that combines the different methods of Sigma and 
HarmonicSums. 
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3.1. All N- Results for 3 Loop Ladder 
Graphs 

Consider, e.g., the diagram containg three mas- 
sive fermion propagators, 



(13) 



From its Feynman integral representation around 
1000 sums have been produced automatically [2"T] 
whose combination leads to an alternative repre- 
sentation in the Mellin variable N. Now the main 
task is to simplify these sums further. One of the 
simplest sums is ([T]) which simplifies with Sigma 
to (IT^I) . A more typical sum (among the 1000 
versions) is 



W V V (-l) s+r 5i(r + 2) 
2^2^ 2^ (N - r)(r + l)(r + 2) 

j=0 g=l r=Q ct=0 v 



mi 



-j + N + s-2 



)(N-3)\(s-iy.a\ 



(-j+N+a+l)(-j+N+<r+2)(-j+N+s+a) ! ' 

Following the same mechanism as described in 
Subsection 12.21 Sigma simplifies this quadruple 
sum to indefinite nested sums. Finally, using the 
HarmonicSums package we transform these sums 
further to an expression involving 145 S-sums 
where the most complicated instances are 

S 2 ,i,i,i(-1,2,~ -1;AT), £2,1,1,1(1, i 1,2; N). 

As it turns out, all 1000 sums can be repre- 
sented in terms of these generalized harmonic 
sums. Combining these expressions leads to a 
rather big expression for diagram (|13j) containing 
in total 533 S-sums. 

At this point we remark that various rela- 
tions among the found S-sums occur. Now we 
could use Sigma again to find a sum represen- 
tation where all the occurring sums are alge- 
braically independent |14j . Since this would take 
quite a while, we use another important feature 
of HarmonicSums following ideas from [18] : it 
contains, e.g., efficient routines and for speed 
up also precomputed tables that produce repre- 
sentations of harmonic sums ([2]) (up to weight 
8, i.e., |mi| + ••• + |m^| < 8) and representa- 
tions of S'-sums (up to weight 6 where the 



Xi are in certain relations, like for instance xi G 
{-2,-1,-1/2,1/2,1,2}) that are algebraically 
independent. Using this feature, all S'-sums can- 
cel, and we end up (within seconds) at a result for 
diagram (|13[) where the following sums remain 

S_ 4 (AT) , S*_ 3 (N) , S-2 (N) , Si(N) , S 2 (N) , S3 (N) , 
Si(N), S_ 3 ,i(JV), S- 2 ,i(A), S 2 ,_ 2 (A), S 2 ,i(A0, 
S 3tl (N),S- 2 , 1 , 1 (N),S 2 ,i,i(N). 

3.2. 3-Loop All N-Results for the N f Con- 
tributions 

E.g., for the diagram with an outer massive and 
an inner massless fermion line, 




(14) 

768 multi-sums have been set up which have been 
simplified to indefinite nested product-sums. One 
of the simple examples is 



N-2 

E 



j(j + l)(j + 2)(A - j)(J - l)! 2 (-j + N - l)! 2 
—j + N — 1 

_ (-N 3 - 5N 2 - AN + 6) (AH) 2 
= (N — 1) 2 N 2 

| 3(A!) 2 (A 3 + 6A 2 + llA + 6) ^ ( 2i ) 

2-^ 



(N - 1)A(2A + 1)( 



2N\ 
N I 



Combining the 768 sums, the final expression 
for (|14[) requires 9 MB of memory and is given 
in terms of 703 indefinite nested product-sums 
which are either harmonic sums, S-sums or of 
the form as given in (J4j . To this end, Sigma 
finds a sum representation where all sums are al- 
gebraically independent against each other. In 
this case, all the unexpected sums cancel, and we 
get the compact form 

_ 20S(l,iV) 4 32(6iV 3 +61iV 2 -21iV+24)S 3 
27(]V+l)(]V+2) H 81Ar a (JV+l)(JV+2) 

16(48JV 5 +746JV 4 +2697iV 3 +2746AT 2 + 1 104JV+240) S 2 
+ 243N-HN+DHW ( 2647V? + mmN " + 21591N5 

+ 52844/V 4 + 74856/V 3 + 66812/V 2 + 305767V + 264o) Si 
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+ 1212357V 4 + 1902357V 3 + 1507587V 2 + 469647V + 2904 J 

/ 40Sf 32(6JV 3 +61iV 2 -21iV+24)Si 

~*~ V 9(iV+l)(]V+2) ^ 27JV 2 (JV+l)(JV+2) 

16(l24Ar 5 +198Af 4 -2387JV 3 -6162JV 2 -3632JV-480) 



81iV 2 (;V+l) 2 (iV+2) 2 °" ^ ^ 2 

32(9JV 3 -623iV 2 +894JV+276) I6OS1 



81]V 2 (jV+l)(jV+2) 



27(JV+l)(JV+2) 



S3 



8(56AT 2 +169iV+112)S4 4(48JV 2 +101JV+96) Sf 
9JV(iV+l)(]V+2) 9JV(AT+l)(iV+2) 

/ 645! _ 128(iV 3 +9iV 2 -10iV-6) \ 

+ I 3(JV+l)(JV+2) 9JV^(JV+l)(JV+2) j 



64S 3i i 64(3iV 2 +7iV+6) 
3(iV+l)(iV+2) 3JV(iV+l)(]V+2) 



"^2,1, 1 



A / SSi(iV) 2 , 16(3iV 3 -iV 2 +30iV+12)Si 
+ S2 ^ 3 (jv+l)(jV+2) 9JV^(JV+l)(iV+2) 



16(3iV 3 +2JV 2 + 17]V+6) 8(4JV 2 +9iV+8)S2 \ 
9N(N+l)'- ! (N+2) 3iV(JV+l)(JV+2) J + 



4485i 



^ \9(N+l)(N+2) 9(iV+l)(]V+2), 

in terms of the algebraically independent set of 
harmonic sums 

Si (TV) , S 2 (TV) , S 3 (V) , 5 4 (V) ,S 2A {N), 
#s,i(iV),S2,i,i(7V). 

4. Conclusion 

We demonstrated how Sigma assists in the 
task to transform multi-sums to indefinite nested 
product-sum expressions by using the summation 
paradigms of telescoping, creative telescoping and 
recurrence solving in the setting of ITS* -difference 
fields, ft turns out that this mechanism works for 
all instances of our performed 2-loop and 3-loop 
calculations. Since these product-sum represen- 
tations extend substantially the class of harmonic 
sums and TS-sums, these algorithms will be use- 
ful for further calculations that might leave the 
solution space of harmonic sums and S'-sums. 

To perform these computations automati- 
cally, a new package called EvaluateMultiSums 
has been developed by the fourth author. 
During these computations, the package 
HarmonicSums [T7] based on ideas of |1I2I6I18[ 
119120) plays an important role to find harmonic 
sum and S-sum representations in optimal form. 
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